Role of negative-energy states and Breit interaction in calculation of atomic 

parity-nonconserving amplitudes. 
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It is demonstrated that Breit and negative-energy state 
contributions reduce the 2.5a deviation [S.C. Bennett and 
C.E. Wieman, Phys. Rev. Lett. 82, 2484 (1999)] in the value 
of the weak charge of 133 Cs from the Standard Model predic- 
tion to 1.7a. The corrections are obtained in the relativis- 
tic many-body perturbation theory by combining all-order 
Coulomb and second-order Breit contributions. The correc- 
tions to parity-nonconserving amplitudes amount to 0.6% in 
133 Cs and 1.1% in 223 Fr. The relevant magnetic-dipole hy- 
perfine structure constants are modified at the level of 0.3% 
in Cs, and 0.6% in Fr. Electric-dipole matrix elements are 
affected at 0.1% level in Cs and a few 0.1% in Fr. 

PACS: 31.30.Jv, 12.15.Ji, 11.30.Er 



Atomic parity-nonconserving (PNC) experiments com- 
bined with accurate atomic structure calculations pro- 
vide constrains on "new physics" beyond the Standard 
Model of elementary-particle physics. Compared to high- 
energy experiments or low-energy ep scattering exper- 
iments, atomic single-isotope PNC measurements are 
uniquely sensitive to new isovector heavy physics |Q. 
Presently, the PNC effect in atoms was most precisely 
measured by Boulder group in 133 Cs ||. They de- 
termined ratio of PNC amplitude -Epnc to the tensor 
transition polarizability (3 for 7Si/ 2 — 6S%/ 2 transition 
with a precision of 0.35%. In 1999, Bennett and Wie- 
man H accurately measured tensor transition polariz- 
ability /?, and by combining the previous theoretical de- 
terminations of the -Epnc [Hr] with their measurements, 
they have found a value of the weak charge for 133 Cs 
Qw = — 72.06(28) CX pt(34)thcor which differed from the 
prediction § of the Standard Model Q w = -73.20(13) 
by 2.5 standard deviations. They also reevaluated the 
precision of the early 1990s atomic structure calcula- 
tions [[|||, an d argued that the uncertainty of the pre- 
dicted £pnc is 0.4%, rather than previously estimated 
1%. This conclusion has been based on a much bet- 
ter agreement of calculated and recently accurately mea- 
sured electric-dipole amplitudes for the resonant transi- 
tions in alkali-metal atoms. 

In view of the reduced uncertainty, the purpose of this 
Letter is to evaluate contributions from negative-energy 
states (NES) and Breit interaction. It will be demon- 
strated that these contributions correct theoretical -Epnc 
and the resultant value of the weak charge by 0.6% in 
Cs. It is worth noting, that due to the smallness of these 



contributions at what had been believed to be a 1% the- 
oretical error in Cs, the previous calculations have either 
omitted or estimated the contributions from Breit 
interaction only partially Q . The main focus of the pre- 
vious ab initio calculations has been correlation contribu- 
tion from the residual Coulomb interaction (i.e. beyond 
Dirac-Hartree-Fock level). In both calculations impor- 
tant chains of many-body diagrams have been summed 
to all orders of perturbation theory. 

This Letter also reports correction due to NES and 
Breit interaction for £pnc m francium. The interest in 
Fr stems from the fact that analogous PNC amplitude is 
18 times larger in heavier 223 Fr compared to Cs @. The 
measurement of atomic PNC in Fr is pursued by Stony 
Brook group ||. 

The quality of theoretical atomic wave-functions at 
small radii is usually judged by comparing calculated and 
experimental hyperfine-structure magnetic-dipole con- 
stants A, and in the intermediate region by comparing 
electric-dipole matrix elements. It will be demonstrated 
that the corresponding corrections to all-order Coulomb 
values are at the level of a few 0.1%. 

The PNC amplitude of nSi/2 — > n'Si/ 2 transition can 



be represented as a sum over intermediate states mP] 
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The overwhelming contribution from parity- violating in- 
teractions arises from the Hamiltonian 



H 



w 



?W/9nuc(r)75 



(2) 



where Gf is the Fermi constant, 75 is the Dirac ma- 
trix, and p nuc (r) is the nuclear distribution. To be 
consistent with the previous calculations the p n uc{ r ) is 
taken to be a Fermi distribution with the "skin depth" 
a = 2.3/(4 In 3) fm and the cutoff radius c = 5.6743 
fm for 133 Cs as in Ref. f§, and c = 6.671 fm for 223 Fr 
as in Ref. 0. The PNC amplitude is customarily ex- 
pressed in the units of 10~ n z(— Qyj/N), where N is 
the number of neutrons in the nucleus (N = 78 for 
133 Cs and N = 136 for 223 Fr ). Atomic units are 
used throughout the Letter. The results of the calcu- 
lations for 133 Cs are £ PNC = -0.905 x lO^Hi-Qw/N), 
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Ref. f§ and £ PN c = -0.908 x l(r u «(-Q w /iV), Ref. §. 
The former value includes a partial Breit contribution 
+0.002 x 10- 11 i(-Q w /N), and the latter does not. Both 
calculations are in a very close agreement if the Breit con- 
tribution is added to the value of Novosibirsk group. The 
reference many-body Coulomb value 



-^PNC 



-0.9075 x 10 _11 i(-Qw/JV) 



(3) 



is determined as an average of the two, with Breit con- 
tribution removed from the value of Notre Dame group. 
For 223 Fr E PNC = 15.9 x lQ- n i(-Q w /N), Ref - 0> this 
value does not include Breit interaction. 

Ab initio relativistic many-body calculations of wave- 
functions, like coupled-cluster type calculations of j|J5), 
to avoid the "continuum dissolution problem" |l0| , start 
from the no-pair Hamiltonian derived from QED [fll|| . 
The no-pair Hamiltonian excludes virtual electron- 
positron pairs from the resulting correlated wave- 
function. If the no-pair wave-functions are further used 
to obtain many-body matrix elements, the negative- 
energy state (NES) contribution is missing already in 
the second order. Recently, it has been shown that 
the magnetic-dipole transition amplitude both in He-like 
ions [[L2| and alkali-metal atoms |13|] can be strongly af- 
fected by the NES correction. The enhancement mecha- 
nism is due to vanishingly small lowest-order values and 
also due to mixing of large and small components of a 
Dirac wavefunction by magnetic-dipole operator. For the 



NES {E < 



the meaning of large and small com- 



ponents is reversed, i.e., small component is much larger 
than large component. The mixing of large positive- 
energy component with small component of NES results 
in much larger one-particle matrix elements, than in the 
no-pair case. The 2mc 2 energy denominators lessen the 
effect, but, for example, the Rb 55*i/ 2 — 65*1/2 magnetic- 
dipole rate is reduced by a factor of 8 from the no-pair 
value by the inclusion of NES |ll|. The inclusion of 
Breit interaction also becomes important and the size 
of the correction is comparable to the Coulomb contribu- 
tion. Just as in the case of magnetic-dipole operator, the 
Dirac matrix 75 in the weak Hamiltonian Eq. (Q) mixes 
large and small components of wavefunctions. Similar 
mixing occurs in the matrix element describing interac- 
tion of an electron with nuclear magnetic moment (hy- 
perfine structure constant A). As demonstrated below, 
the relative effect for these operators is not as strong 
as in the magnetic-dipole transition case, since the low- 
est order matrix elements are nonzero in the nonrela- 
tivistic limit, but is still important. It is worth noting 
that the problem of NES does not appear explicitly in 
the Green's function |5j or mixed-parity j|] approaches; 
however the correction due to Breit interaction still has 
to be addressed. The NES Coulomb corrections have to 
be taken into account explicitly in the "sum-over-states" 
method Q], employing all-order many-body values ob- 
tained with the no-pair Hamiltonian. 

The analysis is based on the Viv-i Dirac-Hartree-Fock 



potential wave-functions, with the valence and virtual or- 
bitals to calculated in the "frozen" potential of core or- 
bitals a. The second-order correction to a matrix element 
of one-particle operator Z between two valence states w 
and v is represented as 
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This expression takes into account the residual (two- 
body) Coulomb, QijkU an d two-body bijki and one-body 
bij = ^2 a biaja Breit interaction. Static form of Breit 
interaction is used in this work. The tilde denotes an- 
tisymmetric combination b^ki = bijki ~ bijik- Subscript 
i ranges over both core and excited states. Note that 
summation over states i and m includes negative-energy 
states. The NES correction to PNC amplitudes arises in 
two circumstances, directly from the sum in Eq. ([!]) and 
in the values of electric-dipole and weak interaction ma- 
trix elements. If the length-gauge of the electric-dipole 
operator is used, the direct contribution of NES in the 
amplitude Eq. (||) is a factor of 10~ 13 smaller than the 
total amplitude, and will be disregarded in the following. 
The numerical summations are done using 100 positive- 
and 100 negative-energy wavefunctions in a B-spline rep- 
resentation 0] obtained in a cavity with a radius of 75 
a.u. 

The breakdown of second-order corrections to matrix 
elements of weak interaction for Cs is given in Table @. 
The all-order values from Ref. Q are also listed in the ta- 
ble to fix the relative phase of the contributions. The ma- 
trix elements are each modified at 0.6-0.7% level. Most 
of the correction arises from positive-energy Breit con- 
tribution, negative-energy states contribute at a smaller 
but comparable level. The contributions from NES due 
to one-body and two-body Breit interaction are almost 
equal ( B W w B { ^) and, in addition, B { + ] w 25 ( j 2) . The 
same relations hold also in francium. The corrections to 
the relevant length-form matrix elements are overwhelm- 
ingly due to the one-body Breit interaction, and are at 
0.1% level. For example, the all-order reduced matrix 
element (6S 1/2 \\D\\6P 1/2 ) = 4.478 Ref. g is increased 
by 0.005, bringing the total 4.483 into an excellent agree- 
ment with experimental value |15| 4.4890(65). Generally, 
the corrections reduce absolute values of the weak inter- 
action matrix elements, and increase absolute values of 
the dipole matrix elements, therefore, their net contribu- 
tions to -Epnc have an opposite sign. Matrix elements 
of weak interaction are affected more strongly, because 
of the sampling of wave- function in the nucleus, where 
relativity is important. 

As demonstrated in Ref. [Q , the four lowest-energy va- 
lence mPi/2 states contribute 98% of the sum in Eq. 
and for the purposes of this work, limiting the sums to 
only these states is sufficient. The corrections to -Epnc 



(4) 
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are calculated first by replacing the weak interaction ma- 
trix elements with the relevant second-order contribu- 
tions and at the same time using all-order dipole matrix 
elements, and second by taking all-order Hyy matrix ele- 
ments, and replacing D with the appropriate correction. 
In both cases the experimental energies are used in the 
denominators. The needed all-order matrix elements for 
Cs are tabulated in Ref . || . The summary of corrections 
to -Epnc is presented in Table ||. The modifications in 
the weak interaction matrix element provide a dominant 
correction. The contribution due to NES in the Coulomb 
part is insignificant, and is already effectively included in 
the reference many-body Coulomb value Ep NC , Eq. (||). 
The reference value -Ep NC is modified by the Breit con- 
tributions by 0.6%, almost two times larger than the un- 
certainty in the Boulder experiment The modified 
value is 



£PN + C B ( 133 CS) 



-0.902(36) x 1Q- L H(-Q W /N) 



A 0.4% uncertainty had been assigned to the above re- 
sult following analysis of Bennett and Wieman jj) . When 



-^pnc * s combined with the experimental values of tran- 
sition polarizability (3 || and Ep^c/fl 0, one obtains 
for the weak charge 

Qw( 133 Cs) = -72.42(28) cxp t(34)thcor. 

This value differs from the prediction of the Standard 
Model Qw = —73.20(13) by 1.7a, versus 2.5a discussed 
in Ref. || , where a is calculated by taking uncertainties 
in quadrature. 

The only previous calculation of Breit contribution to 
PNC amplitude in Cs has been performed by the Notre 
Dame group ||, using mixed-parity Dirac-Hartree-Fock 
formalism. The one-body Breit interaction has been in- 
cluded on equal footing with the DHF potential, but the 
linearized modification to one-body Breit potential due 
to -ffw (Vpnc-hfb in notation of Ref. [Q) has been omit- 
ted. It is straightforward to demonstrate that because 
of this omission, the comparable contribution from two- 
body part of the Breit interaction has been disregarded. 
In units of 10 _11 «(— Qw/N), the result of the present 
calculation for one-body Breit contribution is 0.003 ver- 
sus 0.002 in Ref. [0. Such disagreement is most prob- 
ably caused by different types of correlation contribu- 
tion included in the two approaches. Treating one-body 
Breit together with the DHF potential effectively sums 
the many-body contributions from one-body Breit inter- 
action to all orders, and presents the advantage of the 
scheme employed in Ref. g. However, the dipole matrix 
elements and energies in the sum Eq. (Q) are effectively 
included at the DHF level in the formulation of Ref. [Q , 
in contrast to high-precision all-order values employed in 
the present work. The difference between the two val- 
ues can be considered as a theoretical uncertainty in the 
value of the Breit correction. Clearly more work needs 
to be done to resolve the discrepancy. The accuracy of 
the present analysis can be improved if the one-body Breit 



interaction is embodied in DHF equations, and the many- 
body formulation starts from the resulting basis. How- 
ever, to improve present second-order treatment of the 
two-body part of the Breit interaction, higher orders of 
perturbation theory have to be considered. Apparently 
the most important contribution would arise from terms 
linearized in the Breit interaction, i.e. diagrams contain- 
ing one matrix elements of the Breit interaction and the 
rest of the residual Coulomb interaction. 

The Breit and NES corrections to PNC amplitude 
in heavier Fr are more pronounced. The 223 Fr PNC 
amplitude 15.9 xl0~ n i(-Q w /N) from Ref. @ is re- 
duced by 1.1%. Using all-order dipole matrix elements 
from Ref. the following corrections due to modifica- 
tions in the hw are found (in units of 10~ n i(— Q-yj/N)): 



one-body Breit = —0.131, two-body Breit B± 
—0.053, and the C_ correction, implicitly included in 
Ref. [Q, is —0.003. As in the case of Cs, the all-order 
no-pair Coulomb result || for reduced matrix element 
(JPi/2\\D\\7S\/2) = 4.256 is increased by inclusion of the 
Breit interaction and NES by 0.0011, a 0.3% modifica- 
tion, leading to a much better agreement with experi- 
mental value 4.277(8) |16|. The modification of the £pnc 
due to corrections in the dipole matrix elements is much 
smaller than in the case of /iw- At present there is no 
tabulation of accurate matrix elements of weak interac- 
tion for Fr, and the influence on £pnc due to the Breit 
contribution in dipole matrix elements is estimated from 
average of the modification of individual dipole matrix 
elements 0.2%. The net result decreases the reference 
Coulomb value for 223 Fr § by 0.18 x lQ- u i(-Q w /jV), 
and the corrected value is 



,(2) 



E 



C+B (22Z 
PNC 



Hft) = 15.7 x 10- 1L i(-Q w /N) 



Finally, it is worth discussing Breit and NES contribu- 
tions to hyperfine-structure magnetic-dipole constants A 
for the states involved into PNC calculations. The all- 
order no-pair Coulomb values in the recent work f| have 
been corrected using a similar second-order formulation; 
no details of the calculation have been given. The ex- 
plicit contributions listed in Table III will be useful for 
correcting ab initio many-body Coulomb values. The ta- 
ble presents the contributions for two lowest valence Si/ 2 
and P\/2 states. The calculations are performed using a 
model of uniformly magnetized nucleus with a magneti- 
zation radius R m given in the table. One finds that the 
additional terms reduce values calculated in the no-pair 
Coulomb-correlated approach. For Cs the corrections are 
of order 0.2% for 6S 1/2 , 0.1% for 7S 1/2 , and 0.3% for 
6P%/2 and 7P\/ 2 . The relative contributions to hyperfine 
constants in heavier Fr are larger, accounting for 0.5% of 
the total value for 7S 1/2 , 0.4% for 8Sy 2 , and 0.6% for 
7P 1/2 and 8P 1/2 . 

This work demonstrates that the Breit and NES con- 
tributions are comparable to the remainder of Coulomb 
correlation corrections unaccounted for in modern rela- 
tivistic all-order many-body calculations and hence have 
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to be systematically taken into account. In particu- 
lar, the Breit interaction contributes 0.6% to parity- 
nonconserving amplitudes in Cs and 1.1% in Fr. The 
correction for Cs is almost twice the experimental un- 
certainty and reduces the recently determined || 2.5 a 
deviation in the value of weak charge from the Stan- 
dard Model prediction to 1.7a. Both hypcrfine constants 
and electric-dipole matrix elements are affected at a few 
0.1%. By including NES and Breit correction, the no-pair 
Coulomb all-order dipole matrix elements || for resonant 
transitions are brought into an excellent agreement with 
the accurate experimental values. 

This work was supported by the U.S. Department of 
Energy, Division of Chemical Sciences, Office of Energy 
Research. The author would like to thank W.R. Johnson 
for useful discussions and H.R. Sadeghpour for sugges- 
tions on manuscript. 

TABLE I. Contributions to matrix elements of the weak 
interaction for 133 Cs in units 10 _11 i(-Qw /N). All-order 
no-pair values are from Blundell et al. M. C_ is the cor- 
rection from negative-energy states for the residual Coulomb 
interaction, and and B^ 2 ' are positive/negative energy 
state contributions from one-body and two-body Breit inter- 
action. Notation x[y] — x x 1Q V . 
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2.72[-2] 


-4.74| 
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-1.53 [-4] 


2.96[-5] 


-8.06[-5] 


3.00[- 
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1.54[-3] 
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TABLE III. Contributions to hyperfme-structure con- 
stants in MHz. Column "Expt" lists experimental values, 
where available, and S A gives the total of the contributions 
from negative-energy states and Breit interaction. See the Ta- 
ble | caption for the explanation of other columns. Notation 
x[y] means x x 10 y . 
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m = 5.6748 fm 






6S1/2 


2298.2 


0.11 


-8.14 


0.25 


3.50 


-0. 


35 


-4.64 


7Si/2 


545.90(9) 


0.03 


-1.80 


0.07 


0.96 


-0. 


097 


-0.83 


6P1/2 


291.89(9) 


-6.1 [-4] 


-1.58 


0.25 


0.73 


-0. 


27 


-0.87 


7A/2 


94.35 


-2.2[-4] 


-5.43 


0.09 


0.26 


-0. 


098 


-0.29 
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n = 6. 
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8713.9(8) 


0.07 


-66.7 


-0.08 


19.8 
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54 


-47.4 
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1912.5(1.3) 


0.02 


-12.8 


-0.02 


5.08 


-0. 


14 


-7.88 


7Pl/2 


1142.0(3) 


-5.6[-3] 


-10.8 


1.23 


3.62 
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-6.90 
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362. 91 a 


-2.0[-3] 


-3.61 


0.44 


1.29 


-0. 


34 


-2.23 



a All-order many-body calculations Ref. 



TABLE II. Summary of corrections to PNC amplitude 
in 133 Cs due to Breit interaction and negative-energy states. 
Line SHw lists contributions due to modifications in the weak 
interaction matrix elements, and 8D due to corrections in the 
electric dipole matrix elements. See the Table |j] caption for 
the explanation of columns. The units are 10 _11 i(— Qw/N) , 
and x[y] = x x W y . 
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0.0002 


0.0042 


0.0019 


0.0063 


5D 


-2.6[-10] 


-0.0008 


1.3[-6] 


-0.0008 


Total 


0.0002 


0.0034 


0.0019 


0.0055 
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